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Abstract 

Inaccuracies in the length of members or diameters of joints of large truss reflector 
backup structures may result in unacceptably large surface distortions and member forces. 
However, it may be possible to accurately measure these length or diameter errors. The 
present work suggests that a member and joint placement strategy can markedly reduce 
the distortions and internal member forces. To demonstrate the potential of this approach, 
a member and joint exchange algorithm developed herein is applied to several 102-member 
and 660-member truss reflector structures. It is shown that member and joint exchanges 
can simultaneously reduce the rms surface error and the rms member forces by two orders 
of magnitude. 

* Aerospace Engineer, Structural Mechanics Division, 
f Christopher Kraft Professor of Aerospace and Ocean Engineering. 
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Nomenclature 


Ai = cross-sectional area of i-th member 

Bj = weighting matrix for displacement rms, Eq. (3) 

B s = weighting matrix for member force rms, Eq. (10) 

d = vector of displacements relative to the best-fit surface 

D = reflector diameter 

d rms = displacement rms, Eq. (3) 

d orrna = displacement rms before member exchange 

d rm , = nondimensional surface error parameter, Eq. (28) 

Ei = Young’s modulus of i-th member 

f — vector of external forces 

/ = focal length of the ideal paraboloid 

fi = force simulating i-th member length error, Eq. (7) 

g = d 2 rma gradient vector, Eq. (19) 

H = <C. Hessian matrix, Eq. (17) 

J = objective function combining d rm3 and s rm3 

K = system stiffness matrix 

li — length of i-th member 

P = projection matrix from the best-fit analysis, Eq. (5) 

R = matrix of rigid body modes used in the best-fit analysis 
S = member force influence coefficient matrix, Eq. (15) 
s = vector of internal forces 

s rm , — member-force rms, Eq. (10) 

Sormj = member force rms before member exchange 
s rm , = nondimensional member force parameter, Eq. (29) 

T = coordinate transformation matrix, Eq. (8) 

U = displacement influence coefficient matrix, Eq. (14) 
u = vector of displacements 

w = vector of top surface displacement components normal to the ideal reflector surface 
z i = adjoint load vector associated with i-th component of d 
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a = vector of amplitudes obtained from solution of Eq. (4) 

/ 3 = objective function weighting parameter, Eq. (11) 

8 =minimum acceptable change in the objective function during a member exchange 

A, = adjoint vector associated with i-th component of d 

cr e — standard deviation of the error strain in truss members 

Introduction 

Manufacturing errors in the length of individual members and in the size of joints cause 
two problems in large truss structures. First, many truss structures for space applications, 
such as those supporting antenna reflectors, are required to have an accurate surface shape. 
Errors in member length distort this surface shape (e.g. Ref. 1). Second, because most 
truss structures are structurally redundant, member length errors cause built-in forces in 
the truss members (e.g. Ref. 2). These forces can be a significant fraction of the design 
load for the member and can also adversely affect procedures for assembly of the truss. 

A number of studies have addressed the use of actuators in the structure to reduce 
the surface distortion due to manufacturing errors or other sources (e.g. Refs. 3-5). One 
conclusion of these studies is that many actuators are needed to reduce the surface error 
significantly (e.g. Ref. 5). Surface errors and member forces may also be reduced, without 
resorting to active control, by imposing stricter manufacturing tolerances on members and 
joints. However, this approach increases manufacturing costs and may not be capable of 
achieving the desired surface accuracy for applications demanding extreme precision. This 
paper proposes an alternate approach to dealing with the problem of manufacturing errors 
in truss components. 

Errors in member lengths or joint diameters due to the manufacturing process are 
generally assumed to be random in nature. Previous studies have assumed probability 
distributions for the member length errors and used statistical methods to estimate the 
effects on surface accuracy and member forces (e.g. Refs. 1-3). The output from these 
studies describes the surface error and member forces statistically. 

The present paper, on the other hand, assumes that it is possible and desirable to 
accurately measure the member lengths and joint diameters for a given truss. Then, for a 
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particular arrangement of members and joints, it is possible to predict the specific set of 
surface distortions and member forces of the truss. 

Then the obvious question arises as to which arrangement of members and joints is 
best. This paper develops a method for finding arrangements which reduce the surface 
distortion and member forces. Several numerical studies are performed using tetrahedral 
trusses with 102 and 660 members to demonstrate the effectiveness of this method. In most 
cases the errors are generated as zero-mean random numbers with a normal distribution. 
In one case, the member length errors and joint diameter errors are taken from measure- 
ments of actual truss hardware. In all cases, the member exchange method finds new 
arrangements of members that greatly reduce the surface distortion and built-in member 
forces. 


Tetrahedral Truss Geometry 

The truss structure considered in this study is one that may be used to support a 
precision segmented reflector or antenna surface, and its geometry is shown in Fig. 1. 
The structure has a hexagonal planform and is characterized by the number of rings of 
members in the truss. A typical ring is shown as the shaded region on the five-ring truss 
in Fig. 1. 

Both flat and paraboloidal truss geometries are considered. In the flat case all members 
are designed to have equal nominal lengths. For the paraboloidal geometry, the joint 
locations are obtained by first generating coordinates for a flat truss. Then the change in 
vertical position of the joints on the concave top surface is found directly from the equation 



( 1 ) 


where x and y are the same as for the flat truss, r 2 = x 2 + y 2 , and / is the focal length 
of the paraboloid. The x, y, and z locations for the joints on the convex bottom surface 
are obtained by requiring the lengths of all members in the core to be equal. In the 
resulting reflector geometry, the nominal lengths of of the members in both the top and 
bottom surfaces vary slightly from one another, with some repetition due to the hexagonal 
symmetry. This variation in member length significantly limits the number of members 
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having a given nominal length, and hence the number of possible member exchanges. When 
the members and joints all have their nominal dimensions, the truss has this ideal geometry, 
and it is assumed to be stress-free. Deviations in member lengths or joint diameters from 
this condition cause distortions and residual forces to develop in the structure. In this 
study it is assumed that joint diameter errors can be directly translated into equivalent 
member length errors. 

Calculation of Surface Error and Member Forces 

The truss structure (see Fig. 1) is assumed to support a reflector surface at its upper 
joints. This surface is usually required to be paraboloidal or possibly flat. Due to the 
member length and joint diameter errors the surface will be distorted from this ideal shape. 
A commonly used error measure for the distorted shape is its best-fit root-mean-square 
(rms) deviation from the ideal shape which is defined as follows herein. 

The best-fit surface is defined as the surface for which the mean-square error between 
the distorted truss joints and the ideal surface is minimum. In determining the best-fit 
surface, the ideal surface is allowed to translate along the z axis and rotate about the x and 
y axes relative to the distorted surface defined by the structural joints. Translations along 
the x and y axes are not included in the best-fit analysis because they are undetermined 
for a flat reflector and have only a small effect for shallow, curved reflectors. The resulting 
best-fit displacement, vector d is given as 

d = w — Rq (2) 

where w is the vector of components of joint displacements normal to the undistorted 
reflector surface, R is a matrix of the three rigid-body modes, and a is a vector of am- 
plitudes. The vector a is determined by minimizing the weighted rms value of d given 
as 

d 2 rm , = d T B^d (3) 

The matrix is a positive semidefinite weighting matrix (an identity matrix divided by 
the number of components of d in the examples herein). The vector a which minimizes 
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dims is found by substituting Eq. (2) into Eq. (3) and differentiating with respect to a. to 
get the solution to the normal equations, 

( 4 ) 

so that 

d — [I — R(R T B d R)- 1 R T B (i ]w = Pw (5) 

where P is called the projection matrix. The finite element model of the truss incorporates 
statically determinate rigid body restraints to make the stiffness matrix K nonsingular. 
The finite element analysis, based on small- deformation linear analysis, solves the system 

Ku = f (6) 

where u is the restrained motion due to member length errors and f is the set of equivalent 
member forces that simulate the effect of member length errors. The contribution of the 
i-th member of the truss is a pair of forces colinear with the member of magnitude 

, EiAi 

n = —e< 

where Ei, „4i, and e; are the member Young’s modulus, cross-sectional area, length and 
length error, respectively. 

The vector of upper surface joint displacements normal to the ideal surface is obtained 
by substituting the solution to Eq. (6) into 

d = PTu (8) 

The matrix T extracts from the components of u at each joint in the computational 
coordinate system a displacement normal to the undistorted surface given as 

zu, = tiU (9) 

where W{ is the normal displacement, and u is the vector of displacements in the compu- 
tational coordinate system. The components of the vector t i, the i-th row of T , are all 
zero except for the three direction cosines of the surface normal for the i-th joint at the 



cr = (R T B d R) 1 R T B (i w 
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appropriate locations in the row. Then d rms for the distorted surface can be calculated 
using Eqs. (8) and (3). 

The member force vector, s , is also of interest and can be calculated using the solution 
vector, u . For the member forces a measure of magnitude is 5 rma defined as 

= sT B< s (10) 

where B, is a positive semidefinite matrix (an identity matrix divided by the number of 
members in the examples herein). 

Member Exchange 

Once the members of the truss and its joints are manufactured, an initial configuration 
corresponding to a particular placement of member and joints is selected. The initial 
configuration is defined by a set of member and joint errors with corresponding normal 
displacement and member force measures, d orms and s orma , respectively. A compromise 
objective that combines distortion and member forces is selected as 

2 2 

+d -d)[^A (id 

\“ormj/ \ s orm> J 

The value of (3 can be selected to emphasize either low distortion or low forces or some 
compromise between the two requirements. 

The objective function is then reduced by a binary or triple member exchange algo- 
rithm. The binary member exchange algorithm takes the following form 

1. Select a minimum acceptable objective function improvement fraction 8 (per itera- 
tion), to invoke a member exchange. 

2. Start with the first member. 

3. Calculate the change A J obtained by exchanging the current member with other 
members (starting with the first) until an exchange is found that results in (—A J j J) > 
8 . 

4. Execute the exchange (if any). 

5. If not at the last member, move to the next member and go back to step 3. 
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6. If at the last member, check whether the total number of iterations is exceeded; if yes, 
stop. 

7. If the total number of iterations is not exceeded, reduce 8 and go to step 2. 

Tlie binary exchange strategy is not guaranteed to lead to the minimum. However, 
increasing the number of members involved in a single exchange improves the chances of 
locating the minimum. When all members can be exchanged at each step the minimum is 
guaranteed to be found in a single iteration, but the number of exchanges that is considered 
in that single iteration is the factorial of the number of members. For illustration the 
Appendix presents a simple system where no binary exchange can improve the objective 
function but a triple exchange can reduce it to almost zero. Accordingly, a more expensive 
triple exchange was also used for some examples. The triple exchange algorithm follows the 
same procedure as the binary exchange, except that in step 3 above, the current member 
is matched with all combinations of any two other members. Of course, a triple exchange 
algorithm still does not guarantee that the absolute minimum of J will be found. 

Efficient Reanalysis 

For trusses with several hundred members, the member exchange algorithm requires 
the evaluation of hundreds of thousands to millions of configurations. Therefore, a tech- 
nique for the inexpensive evaluation of the effect of an exchange is required. Because the 
force vector f is linearly proportional to the member length and joint errors, and d is a 
linear function of f , 

d = U;e; 4- U r e r (12) 

where e; and e r are vectors of member length and joint diameters errors, respectively, and 
U( and U r are matrices of influence coefficients. Similarly, the vector of member forces 
can be written as 

s =. S;e; + S r e r (13) 

where Sj and S r are matrices of influence coefficients. The error vectors e; and e r can be 
concatenated into a vector e , and the matrices of influence coefficients U/ and U r and S; 
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and S r similarly concatenated into matrices U and S , respectively. Then Eqs. (12) and 
(13) may be written as 

d — Ue (14) 

and 

s = Se (15) 

Now using Eqs. (3) and (14) the mean-squared error can be written as 

dl ms = e T He (16) 

where 

H = V T B d \J (17) 

so that 

A d 2 rms = g T Ae + Ae T HAe (18) 

where Ae denotes the change in e due to member exchange, and where g, the gradient of 
dims is 

g = 2He (19) 

Note that in spite of its similarity to a quadratic approximation, Eq. (18) is exact, because 
dl ms is a quadratic function. When members i and j are exchanged Eq. (18) gives 

A d r ms ~ (di ~~ 9j){ e j ~~ e i ) + ( e j ~~ e i ) {h-ii + hjj — 2 hij) (20) 

Also, after an exchange the gradient vector g has to be updated by 

Ag = HAe (21) 

In a single iteration Eq. (20) is evaluated a large number of times (up to m{m — 1) times, 
where m is the number of members). The more expensive Eq. (21) is executed only when 
an exchange takes place (up to rn times). Note, however, that for binary exchanges the 
vector Ae has only two nonzero components so that Eq. (21) is inexpensive to evaluate. 
The effect of member exchanges on s rm , are calculated in a similar way. 
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For triple exchanges with members i,j,k permuted to we have 

& d lm, — di{ e k -ei) + gj {ei -ej)+$fc(e,- - e k ) 

-f- hii[c k e^) hjj^Ci e j ) 4" h-kk^-j 

( 22 ) 

T 2hij(e k ^ j ) + 2/i^'^(e{ * ^fe) 

+ 2h ik (ej - e k )(e k - e;) 

Calculation of Influence Coefficient Matrices 

Two techniques are described for calculating the matrices of influence coefficients U 
and S . The first assumes that both the distortions, d , and member forces s are of interest 
and therefore both U and S are required. This technique calculates U and S directly by 
calculating d and s for a unit error in each member. The second technique is based on an 
adjoint variable approach and is substantially cheaper than the direct approach when only 
U is required. 

Direct Method 

From Eq. (12) it can be seen that the i-th column in Uj is just the vector d when a 
unit error is introduced in the i-th member only. This vector can be calculated using Eqs. 
(6), (7) and (8). Each column of the matrix Si is obtained by calculating the member 
forces associated with this displacement solution. The matrices of influence coefficients U r 
and S r can then be formed as linear combinations of the columns in U ; and Si . It should 
be noted that this approach requires the solution of the system equations for each member 
in the truss. For the five- ring truss example, 660 solutions are required. However, the 
stiffness matrix K needs to be factored only once because the effect of the dimensional 
errors affects only the load vector. 
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Adjoint Variable Method 

A substantial savings in computational effort can be achieved when only surface dis- 
tortions are a concern, and therefore only U is required. Since the number of components 
in d is typically much smaller than the number of members in the truss, a dummy load 
or adjoint approach is the preferred method for calculating U . From Eq. (8) the i-th 
component of d may be written as 

=bfd = bfPTu (23) 

where b^ is a Boolean vector with zeroes everywhere except for a one in the i-th position. 
The adjoint approach is employed by noting that the i-th row of TJ i is the vector of 
derivatives of di with respect to the member length errors. Using the adjoint approach 
(e.g. Ref. 6) we apply an adjoint load z i 

*‘ = (^) T = TTpb - < 24 > 

to the truss and solve for the adjoint field A^ 

KAi = z i (25) 

This equation needs to be solved for each component of d . For the five ring case, 91 
solutions are required, as compared to 660 solutions with the direct method. The derivative 
of di with respect to is then found as 

ddi rp di 

d^r l d^ 

The right-hand side of Eq. (26) is the scalar product of the adjoint field by the force 

associated with a unit error in the j-th member which is also equal to the elongation of the 

j-th member due to the adjoint load times the force in that member due to a unit length 
error. This may be written as (see Eq. (7)) 

ddi Ej Aj 

= ^ ~T~~ = 

where e\ tj is the elongation of the j-th member due to A; and s\ t} is the force in the j-th 
member due to A;. Equation (27) implies that the i-th row of U can be calculated by 
applying the load z i of Eq. (24) and calculating the resulting member forces. 
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Numerical Studies 


Several numerical studies have been performed to assess the potential for improve- 
ments in surface accuracy and member forces in tetrahedral trusses. Two- and five-ring 
flat trusses are considered with all members nominally the same length. A five-ring curved 
truss is also studied to assess the effectiveness of member exchanges for this case. Finally 
a study was performed using actual measured properties for a two-ring truss being built 
at NASA-Langley for an assembly experiment. In all cases all members have the same 
cross-sectional area and material properties. 

Except for the assembly-experiment truss, the member length errors are produced 
using a numerical routine for generating random numbers having a normal distribution. 
In these cases it is convenient to generalize the results by defining two nondimensional 
parameters that were also used in Ref. 2. The first is a surface error parameter defined as 

d^T, = (28) 

where D is the reflector diameter and <r e is the standard deviation of the error strain in 
the members. The second is a member force parameter defined as 

^rmi 

&rms — ri a 

EA(T € 

where EA is the member axial stiffness. 

Two-Ring Reflector Example 

A flat, two-ring tetrahedral reflector truss was used as the first example. All 102 
members in the truss have the same nominal length so that all of them can be exchanged. 
The member errors were generated as a set of zero- mean random numbers having a normal 
distribution, and the joint errors were assumed to be zero. After some numerical exper- 
iments, an initial value of the minimum acceptable improvement per exchange, 6 , was 
selected as one percent and this value was reduced by a factor of three at each iteration. 

The iteration history of d rms and s rma is shown in Table 1 for binary exchanges and 
in Table 2 for binary exchanges followed by triple exchanges. The strategy of applying 
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binary exchanges for each iteration followed by triple exchanges for the same iteration 
was used to reduce the computational expense of triple exchanges. It is seen that the use 
of triple exchanges improves the results. For example, for ft 1 the distortion rms is 
reduced by a factor of 146 using binary exchanges and by a factor of 229 using also triple 
exchanges. The corresponding numbers for the force rms (ft = 1) are 237 and 333. Use of 
the compromise objective function (ft — 0.8) results in a reduction of d rrn3 by a factor of 
85 while simultaneously s rm3 is reduced by a factor of 44. 

Convergence is achieved in six or less iterations, and for the first iteration almost all 
members are exchanged while in the last iteration typically only one or two exchanges are 
made (even though the minimum acceptable improvement is greatly reduced). 

Five-Ring Reflector Example 

A flat version of the 660-member five-ring reflector truss shown in Figure 1 was used 
as a second, larger example. All the members in the truss are of the same length so that in 
principle all of them could be exchanged. However, to simulate structures having different 
member lengths, in one case the members were divided into three groups: 240 in the upper 
surface, 195 in the lower surface, and 225 core elements. The first case considered, allowed 
any member to be exchanged with any other, and the second case allowed only members 
within each of the three groups to be exchanged. 

A first set of member length errors was generated using the same procedure as for the 
two-ring case. For this example, only surface accuracy was considered (ft = 1.0). Because 
of the increased number of members only binary exchanges were considered. Again, 8 was 
chosen as 0.01 initially, and was reduced by a factor of three at each iteration. 

The iteration history of d rm3 is shown in Table. 3. Again, convergence is achieved in 
six iterations and d rrn3 is reduced by a factor of 684. Another indication of convergence 
was that for the first iteration almost all members were exchanged, whereas for the last 
iteration only 11 out of 660 members were exchanged. This took place even though in the 
first iteration each exchange required at least a 1 percent improvement in dl ma , while in 
the last iteration the required improvement was 3 6 — 729 times smaller. 
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Table 1. Iteration history for a 2-ring, flat, 102-member truss (random length errors) 
using binary exchanges 


Iteration 

0 = 1.0 

drma ^rms 

0 = 0.8 
dr ms 

&rms 

0 = 
dr ms 

0.0 

&rms 

0 

0.245 

7.55 

0.245 

7.55 

0.245 

7.55 

1 

0.0314 

11.58 

0.0315 

1.94 

0.287 

0.237 

2 

0.410 x 10- 2 

10.12 

0.982 x 10- 2 

0.748 

0.388 

0.0415 

3 

0.199 x 10~ 2 

10.13 

0.520 x 10" 2 

0.424 

0.310 

0.0331 

4 

0.168 x 10“ 2 

10.12 

0.400 x 10~ 2 

0.337 

0.403 

0.0320 

5 

0.168 x 10~ 2 

10.12 

0.262 x 10~ 2 

0.275 

0.408 

0.0318 

6 

0.168 x 10~ 2 

10.12 

0.264 x 10- 2 

0.255 

0.408 

0.0318 

7 

0.168 x 10~ 2 

10.12 

0.263 x 10- 2 

0.255 

0.420 

0.0318 


Table 2. Iteration history for a 2-ring, flat, 102-member truss (random length errors) 
using binary and triple exchanges 


Iteration 

0 = 1.0 

00 

o 

II 


0 = 

0.0 


d rms 

■^rma 

drms 

$rm3 

dr-ms 

^rms 

0 

0.245 

7.55 

0.245 

7.55 

0.245 

7.55 

1 

0.261 x 10- 2 

10.96 

0.484 x lO" 2 

0.298 

0.308 

0.0341 

2 

0.129 x 10- 2 

10.96 

0.293 x 10- 2 

0.174 

0.324 

0.0253 

3 

0.111 x 10- 2 

10.95 

0.287 x 10~ 2 

0.171 

0.312 

0.0253 

4 

0.107 x 10- 2 

10.95 

0.287 x 10- 2 

0.171 

0.278 

0.0239 

5 

0.107 x 10~ 2 

10.95 

0.287 x 10- 2 

0.171 

0.315 

0.0227 

6 

0.107 x 10- 2 

10.95 

0.287 x 10- 2 

0.171 

0.316 

0.0227 


Next, the case of exchanges limited to groups of members was considered. Seven 
iterations were applied to members in the upper surface, followed by seven iterations to 
members in the lower surface, and finally seven iterations to core elements. The results 
shown in Table 4 indicate slower convergence than in the case of unlimited exchanges. 
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Table 3. Iteration history of d rrn , for unlimited member exchanges, five-ring, flat, 660- 
member truss 


Iteration 

(Irma 

0 

0.234 

1 

0.317 x 10- 2 

2 

0.822 x 10“ 3 

3 

0.475 x 10- 3 

4 

0.358 x 10~ 3 

5 

0.343 x 10“ 3 

6 

0.342 x 10“ 3 

7 

0.342 x 10- 3 


Also, the total improvement was a factor of 304. It should be noted that exchanges in core 
elements were most effective. 

Better convergence was obtained when single iterations were performed successively on 
each group. The results, summarized in Table 5, show an improvement factor of 324. The 
effect of changing the order of the groups was also investigated and found to be minimal. 

To check the sensitivity of the results to the set of member errors, another set of 
errors was generated. For this second set the initial value of d TTn , was 0.316. The final 
values were d r m, = 0.216 x 10“ 3 with unlimited member exchanges (a factor of 1460), and 
d rrn3 = 0.815 x 10~ 3 with exchanges limited to each group (a factor of 388 improvement). 
These values are 37 percent lower and 13 percent higher, respectively, than the final values 
for the first set of errors. 

Five-Ring Curved Reflector Example 

This example is similar to the five-ring example above except that the truss is pa- 
raboloidal (as seen in Figure 1) with an f/D of the reflector defined to be 1.0. The main 
consequence of this geometry is that members of many different lengths are now present in 
the truss which limits the opportunities for exchanges. As mentioned above, the geometry 
is generated such that the core member lengths are all equal. Thus, the 225 core elements 
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Table 4. Iteration history of d rms for exchanges limited to groups for a five-ring, flat, 
660-member truss. Iterations completed in one group before the next group was 
started. 


Iteration Upper surface Lower surface Core elements 


0 

0. 

234 


0.428 

X 

10- 

-l 

0.160 

X 

10- 

-l 

1 

0.564 

X 

io- 

-l 

0.277 

X 

io- 

-l 

0.348 

X 

io- 

-2 

2 

0.502 

X 

io- 

-l 

0.223 

X 

io- 

-l 

0.162 

X 

io- 

-2 

3 

0.468 

X 

io- 

-l 

0.197 

X 

io- 

-l 

0.106 

X 

io- 

-2 

4 

0.449 

X 

io- 

-l 

0.184 

X 

io- 

-l 

0.905 

X 

io- 

-3 

5 

0.438 

X 

io- 

-l 

0.176 

X 

io- 

-l 

0.856 

X 

io- 

-3 

6 

0.432 

X 

io- 

-l 

0.168 

X 

io- 

-l 

0.802 

X 

io- 

-3 

7 

0.428 

X 

io- 

-l 

0.160 

X 

io- 

-l 

0.770 

X 

io- 

-3 


Table 5. Iteration history of d rrns for exchanges limited to groups for a five-ring, flat, 
660-member truss. One iteration per group at time. 


Iteration Upper surface 

Lower surface 

Core elements 

0 

0.234 

0.234 

0.234 

1 

0.564 x 10- 1 

0.381 x 10- 1 

0.697 x IO" 2 

2 

0.568 x 10~ 2 

0.431 x IO" 2 

0.213 x 10~ 2 

3 

0.180 x IO" 2 

0.162 x IO" 2 

0.108 x IO" 2 

4 

0.101 x IO" 2 

0.957 x IO” 3 

0.850 x 10“ 3 

5 

0.797 x IO" 3 

0.793 x IO" 3 

0.777 x 10- 3 

6 

0.771 x 10- 3 

0.765 x IO" 3 

0.744 x 10- 3 

7 

0.732 x 10~ 3 

0.724 x 10~ 3 

0.723 x 10~ 3 


plus 48 elements in the faces all have the same nominal lengths and can potentially be 
exchanged with each other. The remaining 387 members are divided into 22 groups, seven 
of which have only 6 members. 

The member length errors were generated as zero-mean random numbers having a nor- 
mal distribution. The exchange algorithm was then run with binary exchanges restricted 
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to members within each group. The objective function was chosen with = 1.0, so that 
only surface distortion was considered. The resulting values of d rrnt for each of the seven 
iterations are shown in Table 6. Despite the constraints on allowed exchanges, the surface 
accuracy, as measured by d rmt , has been improved by a factor of 557. 

Table 6. Iteration history of d rms for the five-ring, curved, 660-member truss. 


Iteration 

drms 

0 

.334 

1 

.164 x 10- 1 

2 

.371 x 1CT 2 

3 

.164 x 10~ 3 

4 

.102 x 10~ 3 

5 

.779 x 10- 3 

6 

.646 x 10- 3 

7 

.600 x 10- 3 


Two-Ring Assembly-Experiment Truss Example 

A precision, two-ring, tetrahedral truss structure (described in reference 7) is be- 
ing built at NASA-Langley for use in an experiment on robotic assembly techniques for 
truss structures. The geometry is nominally flat. The members are constructed from 
2m graphite-epoxy tubes with aluminum end fittings designed to facilitate the automated 
assembly process. During the assembly process, the member is connected to a joint as- 
sembly constructed from a 2.25 inch diameter (nominal), machined aluminum ball with 
attachments to accommodate the member end fittings. 

Manufacturing errors arise from several sources. First, errors in the joint-balls occur in 
the machining process. Measurements of these balls showed them to be highly symmetric 
but with significant differences in the diameter from ball to ball. Figure 2a shows a 
histogram of the measured diameter errors for the 31 balls used in the truss. Errors in the 
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members arise in the fabrication process and in machining operations on the member-to- 
ball attachment fittings. Prior to this study a simple exchange procedure was used to match 
two ball-attachment fittings with each tube with the goal of producing a set of members 
with nearly uniform lengths. Then this combination of tube and two ball-attachment 
fittings was considered as the “member” with an associated length error for use in this 
member exchange study. Figure 2b shows a histogram of the measured distribution of 
length errors for the 102 members used in the truss. 

First, the joint and member errors were arbitrarily assigned in the finite element model 
of the truss. The analysis of this configuration produced a d rmt = .0023 inches and an 
s rma = 25.1 lbs with a maximum member force of 76.7 lbs. From Figures 2a and 2b it can 
be seen that the value of d rm3 is of the same order as the errors in the individual joint 
balls and members. 

The built-in forces in the members are of particular concern in the assembly process 
because they translate directly to the force required by the assembler to insert the member 
into the structure. Accordingly, a value of 0 — .2 was used in the exchange algorithm 
to emphasize the reduction of s rmj over improvement in surface accuracy. Because all 
members in the truss have the same nominal length, any two members are candidates for 
an exchange. The iteration history during the binary exchange process is shown in Table 
7. The final result of d rms = 0.00011 inches and s rms = 0.097 lbs indicates a substantial 
reduction in surface distortion and built-in member forces. These improvements in member 
force and shape accuracy have not yet been duplicated in the experiment, because of 
difficulty in matching analytical and experimental results. 

Concluding Remarks 

Member length or joint diameter errors in truss structures can cause significant shape 
errors or built-in forces in the members. The shape errors are of particular concern when 
the truss structure is being used as a support for a precision electromagnetic reflector. 

A method has been proposed which finds arrangements of members and joints which 
dramatically reduce both surface errors and built-in member forces. The method uses a 
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Table 7. Iteration history for the two-ring, flat, assembly-experiment truss using binary 
exchanges. 


Iteration 

drma (inches) 

Srms 

(lbs) 


Members 

Joints 

Members Joints 

0 

0.231 x lO" 2 

0.231 x 10“ 2 

25.1 

25.1 

1 

0.874 x 10~ 3 

0.645 x lO' 3 

1.47 

1.61 

2 

0.229 x 10~ 3 

0.216 x 10" 3 

0.767 

0.594 

3 

0.154 x 10- 3 

0.148 x 10- 3 

0.235 

0.306 

4 

0.137 x 10- 3 

0.137 x 10“ 3 

0.116 

0.116 

5 

0.115 x 10~ 3 

0.115 x 10~ 3 

0.166 

0.166 

6 

0.112 x lO" 3 

0.112 x lO" 3 

0.094 

0.094 

7 

0.111 x 10‘ 3 

0.111 x 10“ 3 

0.097 

0.097 


heuristic algorithm to select pairs or triples of members for exchanges. Because many eval- 
uations of the effect of potential exchanges on distortions and member forces are required, 
efficient reanalysis techniques based on matrices of influence coefficients of the response 
quantities are used. Both a direct approach and a more computationally efficient adjoint 
approach are used for evaluating the matrices of influence coefficients. 

The member exchange method is demonstrated by application to several two- and five- 
ring, flat and curved tetrahedral truss reflector structures. Sets of member length errors 
were generated as zero-mean random numbers having a normal distribution. For every 
example, the exchange algorithm was able to reduce surface errors and/or member forces by 
at least two orders of magnitude and in some cases three orders of magnitude. The curved 
reflector geometry results in truss members with many different nominal lengths. This 
significantly limits the candidate members for exchange with a given member. However, in 
a five-ring curved reflector example it was still possible to reduce surface distortion by more 
than two orders of magnitude. A two-ring flat truss example using measured data from 
available truss hardware was also considered. In this case also, the exchange algorithm 
was able to substantially reduce computed surface distortions and member forces. 
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Appendix - Benefit of triple exchanges 

The need for triple exchanges may be demonstrated for a simple algebraic problem 
where the matrix H of Eq. (21) is given as 

'2 4-3 

H = 4 10-4 

-3 -4 6.6 

and the error vector e is { — 2 2 7} T so that 

<C. = e T He = 311.4 


It is easy to check that all binary permutations of e increase d rma while the triple permu- 
tation e = { 7 —2 2} r yields d rmt = 0.4. 
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Figure 1.— Parabolic tetrahedral truss reflector structure. 





(b) Member length 
Figure 2 - Concluded. 
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